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Abstract 

We present preliminary improved measurements of the branching fractions of the color- 
suppressed decays B° — > D°h° where h° represents the three light neutral mesons tt°, n and 
ll>. The measurements are based on a data sample of 140 fb _1 collected at the T(4S*) with the 
Belle detector at the KEKB energy-asymmetric e + e~ collider, corresponding to seven times the 
luminosity of the previous Belle measurements. 
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I. INTRODUCTION 



The weak decays B° — > D^°h° [1], where h° represents a light neutral meson, are 
expected to proceed predominantly through internal spectator diagrams, as illustrated in 
Fig. 1. The color matching requirement between the quarks from the virtual W~ and the 
other quark pair results in these decays being "color-suppressed" relative to decays such as 
B° — > D^ + h~ , which proceed through external spectator diagrams. 

Previous measurements of B° decays into D^n , D°r], D°u, and D°p° by the Belle 
collaboration [2, 3], and of B° into D^°tt° by the CLEO collaboration [4], and of B° 
decays into D^n , D^rj, D^°cu, D°r)' by the BaBar collaboration [5] indicate color 
suppressed branching fractions in the approximate range (2-4) x KT 4 . This is substantially 
in excess of theoretical expectations from "naive" factorization models [6-12] in the range 
(0.3-1.7) x 1(T 4 . 

Several approaches to achieving a better theoretical description [8, 9, 13, 14] have been 
developed. They extend upon the factorization approach with consideration of final state in- 
teractions and consequent simultaneous treatment of isospin amplitudes of color-suppressed 
and color-allowed decays. The possibility that similar effects could have dramatic implica- 
tions on the measurement potential of direct CP violation asymmetries in charmless decays, 
together with some degree of discrepancy between the prior Belle [2] and BaBar [5] mea- 
surements provide strong motivation for more precise measurements of the color-suppressed 
decays. 




d + * d d « < d 



FIG. 1: Tree level internal (left) and external (right) spectator diagrams for B — > Dtt decays. 

In this paper we report improved branching fraction measurements of B° decays into 
D°7T , D°rj, D°uj. The measurements are based on a 140 fb™ 1 data sample, which contains 
152 million BB pairs, collected with the Belle detector at the KEKB asymmetric-energy 
e + e~ (3.5 on 8 GeV) collider [15] operating at the T(4S) resonance. This corresponds to 
seven times the luminosity of the previous Belle measurements [2] and almost twice that of 
the earlier BaBar measurements [5]. 

The Belle detector is a large-solid-angle magnetic spectrometer that consists of a three- 
layer silicon vertex detector (SVD), a 50-layer central drift chamber (CDC), an array of 
aerogel threshold Cerenkov counters (ACC), a barrel-like arrangement of time-of-flight scin- 
tillation counters (TOF), and an electromagnetic calorimeter comprised of CsI(Tl) crystals 
(ECL) located inside a super-conducting solenoid coil that provides a 1.5 T magnetic field. 
An iron flux-return located outside of the coil is instrumented to detect K\ mesons and to 
identify muons (KLM). The detector is described in detail elsewhere [16]. 
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II. EVENT SELECTION 



Color-suppressed B° meson decays are reconstructed from candidate D° mesons that are 
combined with light neutral meson candidates h°. The D° mesons are reconstructed in 
three decay modes: K~n + , K^n + n°, and fT~7r + 7r + 7r~ while the light neutral mesons h° 
are reconstructed in the decay modes: 7r° — > 77, rj — > 77, rj — > 7r + 7r~7r° and uj — > 7r + 7r~7r°. 
The invariant masses at each stage of the decay chains are required to be consistent within 
2.5 — 3(7 mass resolution or natural width windows around the nominal masses of the assumed 
particle types. Vertex and mass constrained fits are performed for decays with charged 
products such as the three D° decays and i] — > 7r + 7r~7r°; mass constrained fits are performed 
on the 7r° — > 77 and r\ — * 77 candidates; and vertex constrained fits are performed on 
uj — > 7T + 7r~7r° candidates due to the large natural width of the iv meson. These kinematic 
fits result in improved energy and momenta of the candidate mesons. 

Charged tracks are required to have impact parameters within ±5 cm of the interaction 
point along the positron beam axis and within 1 cm in the transverse plane. Each track is 
identified as a kaon or pion according to a likelihood ratio derived from the responses of the 
TOF and ACC systems and energy loss measurements from the CDC. The likelihood ratio 
is required to exceed 0.6 for kaon candidates. This requirement is 88% efficient for kaons 
with a misidentification rate for pions of 8.5%. 

The photon pairs that constitute tt° candidates are required to have energies greater than 
50 MeV and an invariant mass within a ±3er (er = 5.4MeV/c 2 ) mass window around the 
nominal 7r° mass. 

Candidate rj mesons that decay to 77 are required to have photon energies E 7 greater than 

\E —E 

100 MeV. In addition the energy asymmetry F 71 , F , is required to be less than 0.9. The r? 

Tl ' T2 

candidates are required to have invariant masses within 2.5a mass windows of the nominal 
mass, where a = 10.6MeV/c 2 for the i] — > 77 mode and 3.4MeV/c 2 for the rj — > 7r + 7r~7r° 
mode. If the photons that comprise the r\ — > 77 candidate are found to contribute to 
any n° — > 77 the candidate is excluded. The 7r° decay products of the rj — > it + it~it 
and iv — > 7T + 7r~7r° candidates are required to have CM momentum greater than 200 and 
500 MeV, respectively. The iv candidates are required to have invariant masses within ±31" 
of the nominal mass value, where T is the natural width of 8.9 MeV/c 2 . 

Invariant masses of the D° candidates are required to be within ±2<j of the nominal mass 
where a is 8, 12 and 5 MeV/c 2 for the i^~7r + , K^7T + n°, and K~7t + 7t + tt~ modes respectively. 
The CM momentum of the 7r° in the K~7t + tt° mode is required to be greater than 400 MeV. 



III. B RECONSTRUCTION 

The B° candidates are reconstructed from combinations of D° and h° using the improved 
energy and momenta resulting from the vertex and mass constrained fits. 

Two kinematic variables are used to distinguish signal candidates from backgrounds: the 
beam-energy constrained mass M\, c = yj (E^ cam ) 2 — \ J^Pi I 2 ) an d energy difference AE = 

E* — -Ebcam where E£ CSum is the CM energy, and E* , p* are the CM energy and momenta, 
respectively, which are summed over the D° and h° meson decay candidates. 

The resolution of M bc is approximately 3MeV/c 2 for all modes, dominated by the beam 
energy spread, whereas the AE resolution varies substantially among modes depending 
particularly on the number of 7r° in the final state. Candidates within the broad region 



5 



I AE\ < 0.25 GeV and 5.2 GeV/c 2 < M bc < 5.3 GeV/c 2 are selected for further consideration. 
Where more than one candidate is found in a single event the one with the smaller J2xf/Ni 
is chosen, where xl an d the number of degrees of freedom iVj are obtained from the the 
kinematic fits to the D° and h°. 

A common M^ c signal region of 5.27GeV/c 2 < M^ c < 5.29GeV/c 2 is used for all final 
states. The signal region definitions in AE are mode dependent with |Ai?| < 0.05 GeV for 
uj — > 7T + 7r _ 7r and i] — > 7r + 7r~7r° modes , and \AE\ < 0.08 GeV for ir° — > 77 and i] — > 77 
modes. The event yields and efficiencies presented in the following sections correspond to 
these signal regions. 

IV. CONTINUUM SUPPRESSION 

At energies close to the T(4S) resonance the production cross section of e + e~ — > qq 
(q = u, d, s, c) is approximately three times that of BB production, making continuum 
background suppression essential in all modes. The jet-like nature of the continuum events 
allows event shape variables to discriminate between them and the more spherical BB events. 

The discrimination power of seven event shape variables is combined into a single Fisher 
discriminant [17] whose variables include the angle between the thrust axis of the B candidate 
and the thrust axis of the rest of the event (cos 9t), the sphericity variable, and five modified 
Fox- Wolfram moments [17]. 

Monte Carlo event samples of continuum qq events and signal events for each of the final 
states considered are used to construct probability density functions (PDFs) for the Fisher 
discriminant [17] and cos9 B , where 9b is the angle between the inflight direction and the 
beam direction in the T(AS) rest frame. The products of the PDFs for these two variables 
give signal and continuum likelihoods C s and C q q for each candidate, allowing a selection to 
be applied to the likelihood ratio C = C S /(C S + C q q)- 

Monte Carlo studies of the signal significance N s / ^N s + 1%, where N s and iVj, are signal 
and background yields (using signal branching fractions from previous measurements), as a 
function of a cut on the likelihood ratio C indicate a rather smooth behavior. Although the 
optimum significance is generally in the range 0.6-0.7, a looser cut of C > 0.5 is applied for 
all modes in order to reduce systematic uncertainties. 

For the B° — > D°u mode the polarized nature of the u allows additional discrimination 
against backgrounds to be achieved with an additional requirement of | cos 6^/1 > 0.3, where 
the helicity angle 9hei is defined as the angle between the B flight direction in the u rest frame 
and the vector perpendicular to the ui decay plane in the uj rest frame. 

V. BACKGROUNDS FROM OTHER B DECAYS 

Significant background contributions arise both from color favored decays and from con- 
tributions from other color suppressed decays (crossfeed) B° — > D*°h°. Some backgrounds 
have the same final state as the signal while others mimic signal due to missing or extra 
particles. 

Generic Monte Carlo samples of BB and continuum qq are used to study the background 
contributions in the M bc and AE distributions. The BB event sample excludes the color 
suppressed modes under investigation and associated D*° h° modes. These signal modes for 
each of the decay chains considered and the corresponding D*° h° decays are generated and 
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FIG. 2: Distributions of (a) M hc and (b) AE for B° -» L>°7r°. Points with error bars represent the 
data, the solid line shows the result of the fit and the dotted line represents the signal contribu- 
tion. The crossfeed contributions are represented by the shaded areas. The vertical dashed lines 
represent the signal regions. For (a) the dashed-dotted line shows the continuum-like background 
contribution, with peaking background contributions represented by the small dashed line. For (b) 
the dashed line shows the continuum background contribution, the dashed-dotted lines show the 
B background components contribution. 



reconstructed separately. They are used to estimate the crossfeed contributions to the other 
modes using the branching fractions measured here and by the BaBar collaboration [5]. A 
combined generic Monte Carlo sample weighted according to the effective production cross 
sections and selection efficiencies of qq and BB is also used. 

Figures 2, 3 and 4 show the M bc and AE distributions after application of all selection 
requirements and with the AE signal requirement applied for the M^ c distributions (a) and 
the signal requirement 5.27GeV/c 2 < Mb c < 5.29GeV/c 2 applied for the AE distributions 
(b). The AE signal requirements used are \AE\ < 0.08 GeV for D tt° mode , — 0.08GeV < 
AE < 0.05 GeV for D°r] , and \AE\ < 0.05 GeV for D°u. These regions are indicated by 
vertical dashed lines on the figures. 

The dominant crossfeed contributions to the D°h° decays are found to arise from the 
corresponding D*°h° decays. These contributions peak at same M^ c as the signal but are 
shifted to the lower side in AE. As can be seen from the figures, the crossfeed contribution 
is substantial in the region —0.25 GeV < AE < —0.10 GeV but quite small in the signal 
region. Within the signal region the fraction of crossfeed is less than 10% of the observed 
yield in all cases. 

For the B° — > D°n mode the color allowed B~~ — > D^*'°p~ modes are found to be 
the dominant backgrounds. Non-reconstructed soft it from D*° — > Z}°7r°, photons from 
D*° — > D°7 and 7r~ from p~ — > 7t~tt° produce the same final state as the signal. However 
the missing particles cause a shift in AE with a broad peak centered at approximately 
AE = —0.2 GeV. In order to reduce contributions from this background, events that contain 
B candidates reconstructed as B~ — > D^*'°p~ within the signal region 5.27 GeV < M^ c < 
5.29 GeV and \AE\ < 0.1 GeV are rejected. This requirement reduces the color allowed 
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M bc (GeV/c ) AE (GeV) 

FIG. 3: Distributions of (a) M bc and (b) AE for B° -> D°r] . Points with error bars represent the 
data, the solid line shows the result of the fit and the dotted line represents the signal contribution. 
The crossfeed contributions are represented by the shaded areas. The vertical dashed lines represent 
the signal region. For (a) the dashed-dotted line shows the continuum-like background contribution, 
with peaking background contributions represented by the small dashed line. For (b) the dashed- 
dotted line shows the sum of B background and continuum contributions. 

contribution in the region —0.25 GeV < AE < —0.10 GeV by about 60% it does little 
to reduce contributions in the signal region, but remains useful to facilitate background 
modelling. The Mt, c distribution of these backgrounds is found to contribute at and slightly 
below the the M^ c signal region. 

VI. DATA MODELLING AND SIGNAL EXTRACTION 

Independent unbinned extended maximum likelihood fits to the M^ c and AE distributions 
are performed to obtain the signal yields. The yields from the M^ c fits are used to extract 
the branching fractions, while the yields from the AE fits are used to cross-check the results. 
In most cases the shapes of the signal and background component distributions in M^ c and 
AE are obtained from fits to MC samples. 

The signal models used are the same for all modes, with the M^ c signals modelled with 
a Gaussian function and the AE signals with an empirical formula that accounts for the 
asymmetric calorimeter energy response, known as the Crystal ball line shape [18], added to 
a Gaussian function of the same mean. The signal models are represented in Figures 2- 4 by 
the dotted lines. Fits of the AE distributions to signal Monte Carlo for each final state are 
used to obtain the signal shape parameters; all the Mb c signal shape parameters are allowed 
to float in fits to data. 

The crossfeed contributions in M^ c and AE are studied using a combination of signal 
Monte Carlo samples from all other color suppressed modes, weighted according to the 
branching fractions obtained here or from the Babar measurements [5]. Smoothed his- 
tograms obtained from this combined sample are used as estimates of the crossfeed contri- 
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M bc (GeV/c ) AE (GeV) 

FIG. 4: Distributions of (a) M bc and (b) AE for B° -> . Points with error bars represent the 
data, the solid line shows the result of the fit and the dotted line represents the signal contribu- 
tion. The vertical dashed lines represent the signal region. For (a) the dashed-dotted line shows 
the continuum-like background contribution, with peaking background contributions represented 
by the small dashed line. For (b) the dashed-dotted lines show the sum of B background and 
continuum contributions. 

butions. In the M^ c case the AE signal region requirement results in very small crossfeed 
contributions, which are fixed at the Monte Carlo expectation. For AE there are consider- 
able contributions in the region —0.25 GeV < AE < —0.10 GeV; the normalization of this 
component is allowed to float in the fit. 

Continuum like backgrounds in the Mt, c fits are modelled by an empirical threshold func- 
tion known as the ARGUS function [19]. The small peaking background contributions are 
modelled by a Gaussian of mean and width and normalization obtained by a fit to the BB 
background Monte Carlo M^ c distribution, using an ARGUS function plus a Gaussian. A 
systematic uncertainty of 50% is assigned to the determination of this small background dis- 
tribution. This treatment allows the vast majority of the background to be simply modelled 
with the ARGUS shape leaving a small but less well known peaking background component 
that represents the deviation from the ARGUS shape. 

In fits to data M^ c distributions the ARGUS background function parameters are fixed 
to the values obtained from fits to combined Monte Carlo BB and continuum background 
samples. The signal parameters are free, as are the normalizations of signal and background. 
The small peaking background and crossfeed contributions are fixed at their expected values. 

The AE background distributions in the B° — ► D°r] and B° — > D°cu are modelled using 
smoothed histograms obtained from a combined continuum and generic BB Monte Carlo 
sample. 

For the B° — > D°tt° mode, the shapes of the AE distribution arising from BB and 
continuum background are very different, necessitating separate modelling. The contin- 
uum shape is modelled with a first order polynomial with slope obtained from fits to the 
continuum Monte Carlo sample. The shape of the BB background is modelled with a 
Gaussian function plus a second order polynomial, with parameters determined from a 
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fit to the generic BB Monte Carlo sample. In fits to data the large peak in the region 
—0.25 GeV < AE < —0.10 GeV that arises principally from the color allowed B~ — > D^ Q p~ 
decays is found to be broader than the Monte Carlo expectation, thus all parameters of this 
color allowed Gaussian are allowed to float in the fit. The normalizations of the contributions 
from the remainder of the BB background, the continuum and the signal are also floated in 
the fit, with the small crossfeed fixed as discussed above. 

The results of the M bc and AE fits for the combined modes are presented in Figures 2, 
3 and 4. 



TABLE I: Measured signal region yields and MC estimates of signal region contributions for 
B° — > D°h? for the combined D° subdecay modes and the individual D° subdecay modes. The 
numbers of signal events iV S j ff and continuum like events N^ g obtained from the M^ c fit are listed 
together with their statistical uncertainties. MC estimates of the contributions from peaking B 
background N p ^ and crossfeeds from other color suppressed modes N xrs are listed together with 
their systematic uncertainties. 



Mode 


N sig 


Nbkg 


Npkb 


TV 

1 y xrs 


D°ir° 


637.1 ±33.9 


613.4 ±25.6 


21.1 ± 10.5 


25.9 ±6.5 


D {K7T)7T° 


216.6 ± 17.1 


94.3 ±10.2 


6.1 ±3.0 


6.9 ±1.7 


D°(KTnr°)iT 


184.5 ± 17.6 


191.9 ± 14.0 


4.8 ± 2.4 


7.6 ±1.9 


D (K7nnr)7r° 


260.0 ±21.7 


322.1 ± 18.6 


13.2 ±6.6 


11.5 ±2.9 


IA(77) 


161.3 ± 16.8 


135.6 ± 12.6 


4.5 ±2.2 


10.7 ± 2.7 


Z? (Kvr)r/(77) 


59.8 ±8.9 


21.0 ±4.8 


0.0 ±0.0 


3.7 ±0.9 


D°(Kirir°)r](^) 


51.7 ±8.7 


38.4 ± 6.3 


4.1 ± 2.1 


2.9 ±0.7 


D° (K 7r7T7r) rj (77) 


70.6 ± 10.7 


72.5 ± 9.4 


2.5 ± 1.3 


4.1 ± 1.0 


L>°??(7r 7r7r) 


68.5 ± 10.6 


41.6 ±6.3 


1.2 ±0.6 


2.0 ±0.5 


D (K7t)ri(Tr°Tnt) 


21.4 ±5.3 


7.2 ±2.6 


0.8 ±0.4 


0.5 ±0.1 


D {K7nr°)ri{7T TT7r) 


16.3 ±5.1 


12.6 ±3.5 


1.4 ±0.7 


0.8 ±0.2 


D° (i^7r7T7r)r7(7r 7r7r) 


31.7 ±6.8 


21.6 ±4.5 


0.4 ±0.2 


0.7 ±0.2 


D°7] 


233.1 ± 19.9 


173.9 ± 13.8 


5.5 ±2.7 


12.7 ±3.2 


L>°w(7r°7r7r) 


191.9 ± 17.8 


154.1 ± 12.1 


3.4 ± 1.7 


3.5 ±0.9 


D°(Kir)Lj(ir iTir) 


73.5 ±9.8 


31.8 ±5.7 


3.0 ± 1.5 


0.8 ±0.2 


D°(Kinr )Lu(iT inT) 


54.0 ±8.9 


43.5 ±6.2 


0.0 ±0.0 


1.3 ±0.3 


D°(K7nnr)oo(ir°7nr) 


73.2 ± 11.0 


77.8 ± 8.5 


1.0 ±0.5 


1.4 ±0.3 



VII. BRANCHING FRACTION RESULTS 
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TABLE II: Efficiency from Monte Carlo eMCi correction factor e Dj ^/ MC , and corrected efficiency 
e corr for the combined modes and the individual submodes. This efficiency is for the Id M^ c 
fit sample, with the AE signal region requirement applied. The relative uncertainty is given in 
brackets. 



Mode 
D°(Ktt)tt° 
D°(K 



iriTTT )7T 







£MC £DA/MC 

0.078 ± 0.002(1.9%) 0.938 ± 0.048(5.1%) 

0.177 ± 0.002(1.4%) 0.979 ± 0.045(4.6%) 

0.047 ± 0.001(2.9%) 0.905 ± 0.044(4.9%) 

0.087 ± 0.001(1.5%) 0.979 ± 0.057(5.8%) 



0.075 ± 0.004(5.5%) 
0.173 ±0.008(4.8%) 
0.042 ± 0.002(5.7%) 
0.085 ± 0.005(6.0%) 



D°(KiriT )r](^) 
D° (i^7T7T7r)ry(77) 



0.060 ± 0.001 
0.134 ±0.002 
0.038 ±0.001 
0.064 ±0.001 



2.0%) 
1.6%) 



0.989 ± 
1.032 ±0 
0.954 ±0 
1.032 ±0 



058(5.9%) 
.057(5.5%) 
.054(5.7%) 
.067(6.5%) 



0.061 ±0 
0.139 ±0 
0.036 ± 
0.066 ± 



004(6.2%) 
008(5.7%) 
002(6.2%) 
004(6.7%) 



D 7](7r TTTT) 

D (KTr)7](7r°7nr) 

D°(KTnr°)ri(TT Tnr) 

D (i^7T7r7r)ry(7r 7r7r) 



0.044 ±0.001 
0.096 ± 0.002 
0.028 ±0.001 
0.046 ± 0.001 



2.5%) 
1.8%) 
3.5%) 
2.1%) 



0.940 ± 
0.981 ±0 
0.906 ± 
0.981 ± 



055(5.9%) 
054(5.5%) 
052(5.7%) 
064(6.5%) 



0.042 ± 
0.094 ±0 
0.025 ±0 
0.046 ± 



003(6.4%) 
005(5.8%) 
002(6.7%) 
003(6.8%) 



D°7] 



0.054 ±0.001 



2.1%) 0.971 ±0.057(5.9%) 0.054 ± 0.003(6.3%) 



D°lu(tt°tttt) 
D°(Ktt)lo(tt°tttt) 
D°(KiriT )Lu(iT inT) 
D° (KTT7rir)io (tt^tttt) 



0.028 ±0.001 
0.060 ±0.001 
0.017 ±0.001 
0.032 ± 0.001 



2.9%) 
2.3%) 
3.8%) 
2.5%) 



0.885 ± 
0.924 ±0 
0.854 ±0 
0.924 ±0 



061(6.9%) 
061(6.6%) 
058(6.8%) 
.069(7.4%) 



0.025 ± 
0.055 ±0 
0.015 ±0 
0.029 ± 



002(7.5%) 
004(7.0%) 
001(7.8%) 
002(7i 



TABLE III: Measured branching fractions for the process B° -► D°h° (xlO" 4 ) using separate D° 
subdecay mode samples, as obtained from the M\-, c fit. The branching fractions are listed with 
statistical and systematic uncertainties. 



Mode D°(Ktt) D°(K7nr°) D°(KTnnr) 

D°ir° 2.20 ±0.17 ±0.18 2.11 ± 0.20 ± 0.24 2.75 ± 0.23 ± 0.27 

D°r](^) 1.89 ±0.28 ±0.18 1.73 ± 0.29 ± 0.22 2.39 ± 0.36 ± 0.26 

DVtt^tt) 1.69 ±0.42 ±0.15 1.32 ± 0.41 ± 0.17 2.66 ± 0.57 ± 0.26 

L>°o;(7r 7r7r) 2.55 ± 0.34 ± 0.27 1.95 ± 0.32 ± 0.27 2.45 ± 0.37 ± 0.29 



Results of Mb c and AE fits are consistent; the agreement is typically within 50% of the 
statistical uncertainty. The results from the M^ c fits are found to have a slightly smaller 
total uncertainty in most cases and are used for the final result. Yields are obtained both 
from the individual subdecay mode samples and from samples with the three D° subdecay 
samples combined. The yields from the one dimensional Mb c fits are shown in Table I. 
Both peaking backgrounds and crossfeed contributions in the signal region can be seen to 
contribute substantially less than the extent of the statistical uncertainty on the signal yield. 



11 



The yields obtained are interpreted as branching fractions using the number of analyzed 
BB events, the product of subdecay fractions from PDG [20] corresponding to the decay 
of D°h° into the observed final states and the total selection efficiency The efficiency 
for each mode is first obtained from signal Monte Carlo samples and then corrected by 
comparing data and MC predictions for other processes. For the 7r° reconstruction efficiency 
the correction is obtained from comparisons of i] — > 7r°7r 7r° to r\ — > 77 and to r\ — > 7r + 7r~7r°, 
for data and Monte Carlo. The MC efficiency, correction factor and corrected efficiency are 
presented in Table II. The corrections are obtained from the product of correction factors 
relevant to the final state of each submode. 

The branching fraction results for the individual submodes are shown in Table III. The 
combined submode systematic uncertainties and branching fraction results are shown in 
Tables IV and V, respectively. 

VIII. SYSTEMATIC UNCERTAINTIES 



TABLE IV: Systematic uncertainties of the measured branching fractions for B° — > D°h°, for the 
combined D° submode samples, as estimated for the fit results. 



Category 


D°7T° 


£>y 7 7) 


L>V7r 7T7r) 


D°r] 


D lo(tt tttt) 


Tracking efficiency 


2.6 


2.6 


2.6 


2.6 


2.6 


h° efficiency 


2.7 


4.0 


4.0 


4.0 


5.4 


Kaon efficiency 


f.O 


1.0 


1.0 


1.0 


1.0 


Extra 7T° efficiency 


0.4 


0.4 


0.4 


0.4 


0.4 


Likelihood ratio efficiency 


3.0 


3.0 


3.0 


3.0 


3.0 


MC statistics 


2.3 


2.2 


2.8 


2.4 


3.2 


Peaking background 


1.7 


1.4 


0.9 


1.2 


0.9 


Crossfeed 


1.0 


1.7 


0.7 


1.4 


0.5 


AE resolution 


5.0 


5.0 


5.0 


5.0 


5.0 


Modelling ±ler variations 


2.0 


5.4 


0.2 


10.6 


3.1 


Branching Fractions D°, 7r°, 


rj, to 5.2 


5.2 


5.5 


5.3 


5.2 


Number of BB events 


0.7 


0.7 


0.7 


0.7 


0.7 


Total (%) 


9.5 


11.1 


9.9 


14.4 


11.0 



Systematic uncertainties of the combined modes, estimated for the results based on the M^ c 
fits are summarized in Table IV. Uncertainties on the efficiency correction factors relevant 
to each final state are listed in the Tables together with other uncertainties. For the M bc fit 
the uncertainty from the peaking background, which is fixed at the MC expectation in the 
fit, is obtained by propagating a 50% uncertainty on the normalization of this contribution. 
The crossfeed uncertainty is estimated as 25% of the contribution from this source in the 
signal regions. This accounts for uncertainties on the branching fractions of the crossfeed 
contributions and also differences observed between the floated crossfeed contributions in 
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AE fits and the MC expectation. Uncertainties arising from the background and signal 
modelling used are estimated from the changes in the yields as a result of ±1<7 variations 
on the model parameters. The total uncertainty is obtained regarding uncertainties from 
different sources as uncorrelated. 



TABLE V: Measured branching fractions for the process B° — ► D°h° (xlO~ 4 ) using combined D° 
subdecay mode samples, as obtained from the M\, c fit. The branching fractions are listed with 
statistical and systematic uncertainties. The D°r] result is obtained from a combined sample of 
I>V77) and DV 7r+7r ~ 7r °)- 



Mode 


Branching fraction (x 10 4 ) 


D°ir° 


2.31 ±0.12 ±0.23 


L>V77) 


1.77 ±0.18 ±0.20 


D f](lT lTlT) 


1.89 ±0.29 ±0.20 


D°7] 


1.83 ±0.15 ±0.27 


D°u(ir°inr) 


2.25 ±0.21 ±0.28 



IX. CONCLUSION 

Improved measurements of the branching fractions of the color-suppressed decays B° — > 
D°7i°, D°r] and D°u are presented. The results are consistent with the previous Belle 
measurements. The total uncertainty of the new results is two to three times smaller than 
the previous results, mostly due to the seven times larger data sample. However comparing 
the results with those of BaBar [5] and CLEO [4] indicates an approximately 2a difference, 
with all three branching fractions measured here lower than the previous measurements. 

All the branching fraction results are similar, in the range 1.8-2.4 xl0~ 4 . The large 
values disfavour theoretical predictions based on naive factorization descriptions and indicate 
the need for models including final state interaction effects to satisfactorily describe the 
observations. 
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